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Abstract 

The divergent perturbative expansion of the free-energy density of 
thermal SU(3) gauge theory is resummed into a rapidly convergent 
series using a variational implementation of the method of confor- 
mal mapping of the corresponding Borel series. The resummed result 
differs significantly from non-perturbative lattice simulations and the 
discrepancy is attributed to the presence of a pole on the positive axis 
of the Borel plane. The position of that pole is determined numeri- 
cally and the difference between the lattice data and the resummed 
series is related to a phenomenological bag 'constant'. 



It is generally believed that the collision of heavy-ions at sufficiently high 
energies will lead to the formation of a new phase of matter, the quark- 
gluon plasma, and experiments to produce such a plasma are underway at 
Brookhaven and CERN. As the effective coupling, a, of quantum chromo- 
dynamics (QCD) decreases with an increase in energy, theorists have used 
perturbative methods to study properties of the plasma at high temperature 
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(T). For example, a completely analytical calculation of the free-energy den- 
sity of QCD to order a 5 / 2 was performed a few years ago [1], |2|]. The purely 
gluonic contribution is given by 

where Fq = ~ 8 ^ 5 T is the contribution of non-interacting gluons and ft is the 
renormalization scale in the MS scheme. Unfortunately (||) is an oscillatory, 
non- convergent, series even for a as small as 0.2, which is close to the value 
of physical interest. 

Pade' Approximants (P.A.'s) were used in Ref. to resum the series (p]). 
It was found that the dependence on the arbitrary scale p, was reduced and 
the convergence of the series somewhat improved. However P.A.'s have a 
number of well-known drawbacks. For those and other reasons, in Refs.Q 
the authors abandoned the expansion of the free-energy density with respect 
to the coupling constant and considered instead selective resummations of 
gauge-invariant diagrams. Though the results of |3J compare favorably with 
lattice simulations f5j, calculations beyond leading order are complicated 
and thus it seems that the issue of convergence is left open. In Ref. 0, 
yet another procedure was used to study the free-energy of hot QCD: short 
distance perturbative effects were handled analytically while long-distance 
effects were described 0] by an effective three-dimensional theory and studied 
numerically. 

It appears that a number of separate issues concerning Eq.(P have be- 
come confused in the literature. The first issue is whether the given divergent 
series can be resummed into a convergent series, preferably in a systematic 
and well-motivated way. The second issue is whether such a resummed series 
accurately represents the physical quantity. The third issue is whether one 
obtains any new physical insight in the process. It will be the attempt of 
this paper to shed some light on these and related questions. 

Recall that the divergence of perturbative expansions in quantum field 
theory is is a generic phenomena [0] . Given a series 

N 

S N (\) = l + Y,fn\ n , (2) 
n=l 
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where A is the coupling constant, one expects the coefficients f n to grow as 
n! for large n . It is natural then to introduce the Borel transform 

N f n 

5(.) = l + £^- (3) 

n=l U - 

which has better convergence properties than (0). The series (|2]) may be 
recovered from @ through a Laplace transform 

1 r°° 

S N (\) = - dz e- z/x B{z) . (4) 



A Jo 

In order to proceed non-trivially, one first performs an approximate sum- 
mation of the series (^) so that Eq.(|4]) then gives the resummed version of 
(FJ). Now, suppose that the only singularity of B(z) is at z = —1/p, with 
p real and positive. Then the radius of convergence of the Borel series is 
1/p. In order to perform the integral in (f|), one needs to extend the domain 
of convergence of the Borel series. One way to do this is by the method of 
conformal mapping M, H|. Define 



y/TTpZ - 1 
W(Z) = ; 5 

which maps the Borel plane to a unit circle. The inverse of (|5]) is given by 

Z = ~p-Ji^y (6) 

The idea is to rewrite (|3]) in terms of the variable w. Therefore, using 
(H), z n is expanded to order N in w and substituted into The result is 



where is given by (|^). Equation (0) represents a resummation of the 
original series (0). This technique has been used in determining critical 
exponents in statistical systems where the singularity at z — —1/p is due to 
instantons ]7|, |J, and in QCD at zero temperature where the singularity is 
due to renormalons |§. 

Currently no information is available about the exact location of singular- 
ities in the Borel plane of thermal QCD though undoubtedly there is at least 
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one on the negative semi-axis. Therefore in order to apply the resummation 
(0) to Eq.(0), a new idea is introduced in this paper: It is first assumed that 
the only singularity is at z = —1/p, p > 0, with the value of p determined 
by the condition that it be the position of an extremum of ([?]). That is, p is 
chosen to be a solution of 



' ds N (\, P y 

x dp j 



= 0. (8) 

A=A 

Since S^r(X,p) depends on the coupling A, one first fixes A at some ref- 
erence value A = Ao (say, at the mid-point of the range of interest) in order 
to solve Eq.(|). Fortunately, it turns out that the solution of (§), and hence 
the convergence of Sn, is not very sensitive to the exact value of A . 

Let me illustrate the technique by applying it to two cases where exact 
results are known. Consider first the integral 

If the right-hand-side of (§) is expanded as a power series in A one obtains 
at iV-th order 

N 

I N (X) = 1 + E A n (-l) n n! . (10) 

71=1 

Clearly, from (|9]), the exact Borel transform of this series is B(z) = 1/(1 + z) 
with a singularity at z = —1. Ignoring this information, let us resum the 
divergent series fllCf ) using ([7[) with /„ = (— l) n n\ and values of p determined 
for each N through equation (§) at the reference value A = 0.5. The results 
for (|5J) are as follows: There is no extremum for N — 1. For N = 2, there 
is a minimum at p = 2.65. For N = 3 there is a local maximum at p = 1.5 
and a minimum at p = 5.1. For N = 4 there is a local minimum at p = 1.3, 
a local maximum at p — 2.3 and a global minimum at p = 8.4. Jumping 
ahead to iV = 8, there is a local minimum at p = 1.075, a local maximum at 
p = 1.5 and a global minimum at p — 3.6. 

Thus in general (|8|) has more than one solution for a given iV and Ao- In 
Fig.(l), Eq.(|^) is plotted for iV = 2,3,4 and 8, at the respective minimum. 
Notice how the curves converge rapidly to the exact value given by @. 
Alternatively, one might choose for each N, the value of p (from the multiple 
solutions of Eq.(HD) which seems to be part of a converging sequence. In this 
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case the values, p = 2.65(iV = 2),p = 1.5(JV = 3),p = 1.3(N = 4), and 
p = 1.075(A r = 8), appear to converge to the exact value p = 1. The curves 
are shown in Fig. (2). Clearly the curves in Fig. (2) converge faster to the 
exact value than those of Fig.(l), but unless one is interested in very high 
numerical accuracy, the difference is not significant. For example, at A = 0.5, 
the exact value of @ is 0.722657, while the resummed value for N = 8 is 
given at the global minimum p = 1.075 by S$(X = 0.5, p = 1.075) = 0.722652, 
and at the local minimum p = 3.6 by Ss{X = 0.5, p = 3.6) = 0.722524. 

The main points illustrated by this example, which seem to be common 
to the other cases studied, are, (i) the rapid convergence of the resummed 
series represented by (^) compared to the original wildly oscillating series 
(|T(]|), (ii) the relative insensitivity of the convergence and numerical value of 
the resummed series to the particular extremum chosen among the possible 
multiple solutions of (§) (for a given iV and Ao), even if the chosen value 
of p is quite different from the exact value, (iii) the relative insensitivity of 
Eq.(||), and hence the convergence of (0) to the precise value of A . 

For another example, consider thermal O(M) X 2 ^ field theory in the limit 
M — > oo. The exact free-energy density in this case has been determined in 
Ref.||10||. In Eq.(5.8) of that paperQ the perturbative expansion, in A, of the 



free-energy density is also given up to A 6 . Defining S = (F(T) — F '(0))/ 'F i( i ea i 
and choosing JI = T for simplicity, the values of f n for 2 < n < 6 can be 
read off from Eq.(5.8) of Ref.|l0| and (§) solved at some reference value, 
say Ao = 4. The solutions of (|8|) in this case are: N = 3,p = 0.1 (min); 
N = 4,p = 0.05(local max), p = 0.2(min); N = 5,p = 0.025 (local min), 
p = 0.1 (local max), p = 0.3(global min); iV = 6,p = 0.1 (local max), 
p = 0.45 (global min). As in the first example, the convergence of the 
resummed series is found to be rapid even for large coupling, in contrast to 
the oscillatory behaviour of the ordinary pertubation expansion observed in 
HT0| . Fig.(3) shows the curves for Sn(X),3 < N < 6, for the value p = 0.1, 
which seems to be the value p converges to as N increases. The exact value 
of the free-energy density at A = 8 taken from Fig. (6) of Ref.[|l(| is about 
0.875. By comparison the resummed value predicted here is given at sixth 
order by Sq(X = 8,p = 0.1) = 0.889. On the other hand, if one evaluates 
S*6 not at p — 0.1 but rather at its global minimum p = 0.45, one gets 
S 6 (\ = 8,p= 0.45) = 0.849, a difference of less than 5%. 

1 Note that the definition of the coupling constant used here is different from that in 
Ref.@. 



5 



Actually, no information is available about the singularities in the Borel 



plane for the 0(M — > oo) scalar field studied in |H|. The good agreement 
of the results obtained here with the exact results of Ref.fpHfl leads one to 
conjecture that for the free-energy of this model, the singularity closest to 
the origin in the Borel plane might be near p = 0.1, that is, z = —10. 

Note that although the coupling constant in the scalar field theory model 
is A 2 , the perturbative expansion of the free-energy density contains the odd 
powers A 3 and A 5 which is typical of thermal field theories with massless 
particles (or at very high temperatures) |], [TT| . Physically these are due 
to collective effects such as Debye screening and it is sometimes suggested 
in the literature that these terms should be treated on a different footing. 
However, as the analysis above shows, from a mathematical point of view 
these odd powers are no different from the other terms in the expansion of 
the free-energy density and can be resummed as part of a single series. 

Finally, the resummation technique of Eqs.([7|||) is applied to the free- 
energy density of SU(3) gauge theory given in Eq. ([!]). As in Refs.j|], I 
replace (-)a by the approximate two-loop running coupling constant defined 
by 

AM)= , 2 fi- »^'») (ii) 



where L(c,x) = ln((2.287rcx) 2 ), c = ft/2nT and x = T/T c , with T c ~ 
270MeV the critical temperature which separates the low and high tem- 
perature phases |J. Furthermore, as in Ref. 0, I have absorbed the In (a) 
term which appears at three-loop order into the coefficient of the a 2 term in 

©■ 

Fixing first the reference values cq = 1 , xq = 3 (which fixes the reference 
value of Ao), the results of (H) are: N = 2, no extremum; N = 3, p = 3.2(min); 
N = 4, p = 7.6(min); N = 5,p = 13.1(min). Since c and x appear in ([!]) 
and ( |TTD only logarithmically, changing these values in the range of say, 
0.5 <c<2, 2 < x < 5, has almost no impact on the solution of Eq.(H) and 
hence on the optimal values of p. 

The curves for the resummed series are plotted in Fig. (4) for c = 1, that 
is at the renormalization scale p, = 2ttT. Again the rapid and monotonic 
convergence is manifest, the result approaching the ideal gas value even at 
moderate temperatures ~ 2T C . One can estimate the effect of the unknown 
higher order, A 6 , contribution. It turns out to be negligible [|12|]. Therefore 



one feels confident that the N = 5 curve in Fig. (4) is numerically close to 
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the (unknown) exact sum of the perturbation series. In Fig. (5) the curve 
for S 5 (x,c,p = 13.1) is plotted for three values of c to indicate its mild 
dependence (less than 1%) to the arbitrary renormalization scale /2. 

Lattice results for the free-energy density of pure SU(3) theory are shown 
in Fig(6). The lattice community has indicated that their errors are under 
control (less than 5%). In that case, I am left with the task of explaining the 
significant difference (e.g. ~ 15% at T = 3T C ) between the best analytically 
resummed result represented by the N = 5 curve in Fig. (4) and the lattice 
data. At zero temperature, it is known that non-abelian gauge theories are 
not Borel summable pi. That is, B(z) contains singularities for positive z, 
rendering the integral in (Q) ambiguous. The situation is not expected to be 
different at non-zero temperature. Usually 0, the presence of such singular- 
ities is taken to indicate the existence of non-perturbative corrections. One 
can estimate the ambiguity, 5S, and hence the non-perturbative correction, 
as the residue of the integrand in (f|) at the location of the singularity 0. If 
the singularity of B(z), on the positive semi-axis, closest to the origin is a 
pole at z = q, then from (|) 

5S = ^ e' q/x (12) 
A 

where A is a constant. Assuming that the difference between the lattice data 
and Fig. (5), is due to (0), the constants A and q can be determined by 
rewriting fli~2]) as 

\n(\5S) = \n(A) - q/X (13) 

and using for 5S the difference between Fig. (6) and the median curve in 
Fig. (5) (i.e. c = 1). Fig. (7) shows the left-hand-side of ( |i~3|) plotted against 
1/A. This gives A = e 8 ' 7 and q = 2.62, and so, with A(x) = A(c = l,x), 

Su* = S pert - A- e 8 - 7 ' 2 ^ , (14) 
A{x) 

where Si att represents the lattice data for the free-energy, and S pert the Borel 
resummed perturbative result, both normalized with respect to the ideal gas 
value. 

It is extremely reassuring that both the sign and magnitude of q deter- 
mined in this way are self-consistent with the assumptions made. In partic- 
ular, the singularity at z = q = 2.62 is more than 30 times away from the 
origin than the singularity at z = — 1/p = —1/13.1 and justifies a posteriori 
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the resummation procedure (0) which considered only the nearest singular- 
ity. Furthermore since S pert is extremely close to the ideal gas value, Eq.(|14D 
may be interpreted as a generalisation of phenomenological equations of state 
TTj| for the free-energy where the second term on the right-hand-side of ( |Tl| ) 



is called a 'bag constant'. In our case the 'constant' is really temperature 
dependent and represents a non-perturbative contribution to the free-energy 
that vanishes at infinite temperature. Note that this interpretation of the 
second term on the right-hand-side of (|14D is consistent with the usual one 
only because the resummed perturbative result lies above the lattice data. 

Let me now summarize the main results of this paper. A new procedure, 
a variational version of the well-known conformal- mapping of Borel series, 
was introduced and its practicability illustrated. It was shown that the badly 
divergent series for the free-energy density of thermal SU(3) gauge theory ([I]) 
could be resummed in a systematic and relatively well-motivated way into 
a rapidly convergent series. However the final result differed significantly 
from non-perturbative lattice data, suggesting that the discrepancy is due 
to the non-Borel-summability of the theory. Numerically, the difference was 
parametrised in terms of two constants (see (|T4l) ), and it is suggested that the 
ambiguity in the Borel integral is the bag 'constant' used in phenomenological 
models for the free-energy density. 

In physical terms, the situation may be described as follows. If the Borel 
resummed perturbative expansion had agreed with the lattice data, then 
one would have argued that the high-temperature phase of SU(3) theory is 
appropriately described by weakly coupled gluons. However, if the deviation 
of lattice data from the resummed and convergent perturbation expansion 
found here is taken at face value, then one is led to conclude that even at 
temperatures as high as three times T c (~ 700MeV), the phase of thermal 
SU(3) is not accurately described by weakly coupled gluons but rather by 
more complicated structures which are responsible for the non-perturbative 
bag contribution. If one accepts this conclusion, then one must also be open 
to the possibility that what will be produced at CERN and Brookhaven might 
be better characterised as something other than a quark-gluon plasma. For 
some alternative descriptions of the high-temperature phase of QCD see, for 
example, JTB| ] . 



An extension of the analysis presented here to QCD and other thermal 
gauge theories, a further development of the methodology itself and its ap- 
plications to other physical problems, will be presented in an accompanying 



series of papers [12 
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Figure Captions 

Figure 1: Plot of Eq.(0) for the model in Eq.(|g), for N = 2,3,4 and 8 
at the respective global minimum of Eq.(||). The curves move upwards with 
increasing N. The curves for iV = 3, 4 practically coincide, while the curve 
for the exact expression @ is indistinguisable from that for N = 8. 

Figure 2: Plot of Eq.® for the model in Eq.(0), for N = 2,3,4 and 8 
at the respective values of p from Eq.(|8]) which converge to 1. The values of 
the curves at A = 0.5 are, N = 2(0.704), N = 3(0.726), iV = 4(0.7219), N = 
8(0.7227). The curve for the exact expression (||) is hardly distinguisable 
from that for N = 3, 4, 8. 

Figure 3: Plot of Eq.fl?]) for the model in Ref.0 for iV = 3,4,5 and 
iV = 6 at the value p = 0.1. The curves move upwards as iV increases. 



For N = 5, 6 they are practically identical. The exact curve in Ref. [10| lies 



slightly below that for N = 6: At A = 8 the exact curve has the value 0.875. 

Figure 4: The resummed (|7|) free-energy density of hot SU (3) gauge the- 
ory for N = 3,4 and 5, at the renormalization point ]x = 2-kT . The curves 
move upwards as increases. 

Figure 5: The fifth order resummed free-energy density (0) of hot SU(3) 
gauge theory at p = 13.1, for the renormalization scale values fl = 0.5, 1 and 
2. The free-energy density increases with increasing fx. 

Figure 6: Mean lattice results for the free-energy density of hot SU(3) 
gauge theory from Ref. || . Here S[ aU refers to the free-energy divided by the 
free-energy of an ideal gas of gluons. 

Figure 7: Plot of the left-hand-side of Eq.([L3|) against 1/A at several 
points (temperatures). 
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